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We consider properties of zero and near-zero fermionic modes in lattice gluodynamics. The modes are 
known to be sensitive to the topology of the underlying gluonic fields in the quantum vacuum state of the 
gluodynamics. We find evidence that these modes are fine tuned, that is exhibit sensitivity to both physical 
(one can say, hadronic) scale and to the ultraviolet cut off. Namely, the density of the states is in physical units 
while the localization volume of the modes tends to zero in physical units with the lattice spacing tending to 
zero. We discuss briefly possible theoretical implications and also include some general, review-type remarks. 
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1. TWO SCALES OF QCD 

In this note we will consider properties of vacuum 
fluctuations within lattice formulation of Yang Mills 
theories. For simplicity, we concentrate on the case of 
pure gluodynamics, with no dynamical fcrmions. What 
is specific for the lattice formulation (see, e.g., PP) is 
that it is a field theory in Euclidean space-time. The 
action reads: 



(1) 



X 



where G" (x) is the non-Abelian field strength tensor 
and a is the color index, g is the coupling. We will 
consider actually the SU(2) case, a = 1,2,3. Using 
action one generates vacuum field configurations, 
{A^}, where is the gauge potential and performs 
further measurements on these fields. 

As any renormalizable theory, quantum gluodynam- 
ics exhibits two scales, infrared and ultraviolet. More- 
over, the ultraviolet cut off is introduced explicitly, 
through a finite lattice spacing a. The infrared scale, 
Kqcdi on the other hand is emerging dynamically: 
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exp( - b /g 2 (a)) 



(2) 



where bo is a constant, g(a) is the bare coupling con- 
stant normalized at the lattice spacing, Aqcd charac- 
terizes the scale where the running coupling is of order 
unit. As a — > 0, the bare coupling squared, g 2 (a) tends 
to zero as an inverse log of (a ■ Kqcd)- If one changes 
the lattice spacing a and modifies g (a) according to the 
rules of the renormgroup the scale Aqcd does not de- 
pend on a. 

Lattice simulations allow to study directly vacuum 
fluctuations and the both scales, Aqqd and l/a get 



manifested in the vacuum fluctuations. In particular, 
the zero-point fluctuations are sensitive to the ultravi- 
olet cut off. One can measure them, for example, by 
studying the gluon condensate, or the vacuum expecta- 
tion value of the gluonic field strength tensor squared: 
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{1 + E k a k g 2k {a)} (3) 



where a k are coefficients of the perturbative series. The 
matrix element © is divergent as the fourth power 
of the ultraviolet cut off. This is the well known di- 
vergence of the density of the vacuum energy in field 
theory, which arises due to the zero-point fluctuations. 
High gluonic frequencies dominate the matrix clement 
© because of the phase space associated with these 
fluctuations. One can say that zero-point fluctuations 
represent an example of entropy dominated fluctuations. 

Situation looks absolutely different if one consid- 
ers topological charge. The density of the topological 
charge is given by: 



Q(x) 
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where tavpa is t ne totally antisymmetric tensor. One 
can readily show that Q(x) is a full derivative. As a re- 
sult, all perturbative fluctuations do not contribute to 
the topological charge. On the other hand, probability 
to find a non-perturbative fluctuation of size p is sup- 
pressed for small p as exp(— const/ g 2 (p)) where g(p) is 
the running constant. This factor grows fast with grow- 
ing p and is somehow smoothened out at p ~ ^qcd 
where the running of the coupling cannot be calculated 
reliably. As a result, topologically nontrivial fluctua- 
tions have typical size of order and are absolutely 
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negligible on the scale of a. In particular, instantons are 
topologically non-trivial and all the factors mentioned 
above are known explicitly. Instantons represent fluc- 
tuations whose probability is determined primarily by 
action, not entropy (for a review of instantons, in the 
context of the lattice measurements see, e.g., 0). 

A common viewpoint is that confinement is due to 
vacuum fluctuations on the scale Aqcd, i.e. is deter- 
mined by soft, scmiclassical fields. The expectations can 
be confronted with lattice measurements on the fully 
quantum, vacuum state of gluodynamics. And most re- 
cently there has been emerging evidence (Hj that the 
confining fields are actually of a third type, exhibiting 
both infrared and ultraviolet scales. Both action and 
entropy are very large but balance each other almost 
exactly. One can call them fine tuned fluctuations, for 
a review see 0]. 



2. FINE TUNING 

Imagine that a relativistic system has size ro. Then, 
typical momenta should be of order p ~ 1/ro. Re- 
spectively, one expects that typical masses are of order 
m 2 ~ l/ r o- 1^ on the other hand, observed masses are 
much smaller than r^~ 2 , one calls such a case fine tuning. 

The notion of fine tuning has been discussed most 
thoroughly in connection with the Higgs physics. The 
problem here is that to fulfill its role as a part of a renor- 
malizable theory of weak interactions the Higgs particle 
should have mass of order 100 GeV, like intermediate 
bosons of weak interactions. On the other hand, if Hig- 
gses are point-like particles down to scale ro, then ra- 
diative corrections to the Higgs mass are of order 

Am 2 H ~ ar^ 2 , (5) 

where a is the electromagnetic coupling. If ro is, say, of 
order of inverse Planck mass, fine tuning of the radiative 
correction and of the bare mass is required. 

To resolve the puzzle of fine tuning, if it is observed 
experimentally, one usually invokes hidden symmetry. 
In the Higgs case, for example, one of the favorite can- 
didates for such a symmetry now is supersymmctry. 

In case of Yang-Mills theories, one usually does not 
expect to confront the problem of fine tuning. However, 
recently it was found that the field configurations which 
are responsible for the confinement appear fine tuned. 
In particular, there is an ample evidence that so called 
central vortices, for review see, e.g., 0, are responsi- 
ble for the confinement. The central vortices represent 



closed two-dimensional surfaces whose total area scales 
in physical units: 

A tot a 24 Vtot fm- 2 (6) 

where Vtot is the total volume of the lattice. On the 
other hand, the non-Abelian action associated with the 
surfaces is ultraviolet divergent [^j: 

Sua « 0.55 . (7) 

cr 

Combining observations © and (JJJ one concludes that 
the suppression of the fluctuations due to the action Q 
is to be nearly compensated by enhancement due to the 
entropy. That is, the fluctuations are fine tuned. 

The observations above are pure empirical. Theo- 
retically, it is natural to speculate again that there is a 
hidden symmetry which ensures the observed fine tun- 
ing. Moreover, the only symmetry which can come into 
consideration is the conjectured duality between Yang- 
Mills theories and string theories However, such a 
connection is pure speculative at the moment 4 . We 
mention this possibility just to emphasize that further 
studies of the fine tuning in Yang-Mills theories are of 
great interest. In this paper, we address probable man- 
ifestation of the fine tuning in fermionic zero modes. 

3. WITTEN-VENEZIANO, BANKS-CASHER 
RELATIONS 

We will study properties of low-lying modes of the 
Dirac operator, 

£>m7m V*n = ^ni>n , (8) 

where 7^ are the Dirac matrices and the covariant 
derivative is constructed on the gauge potential 
generated as a vacuum field configuration. The Dirac 
equation is solved numerically. Moreover, one generates 
many configurations and studies the properties of the 
modes with low values of A n . Note that both X n and 
the volume occupied by the n-th wave function ip n are 
gauge invariant quantities. 

There exists a rich literature on the low-lying Dirac 
eigenmodes (LDEs). Most commonly, one uses the in- 
stanton model of the vacuum, for review and further 
references see, e.g., [2]. The instantons are localized so- 
lutions with non-trivial topological charge and there are 
exact zero fermionic modes associated with them: 

n + - n_ = Nf ■ Qtop (9) 

where n± are the number of zero modes with positive 
or negative chirality, Nf is the number of quark flavors 
and the topological charge Qtop = 1 for an instanton. 
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In the physical vacuum, the instanton solutions are 
distorted. Indeed, there are neighboring instantons 
(or anti-instantons) and the instanton fields are mod- 
ified because of that. The corresponding fermionic zero 
modes are becoming near-zero modes. 

There are some generic features which are predicted 
to survive the modifications. First, consider exact zero 
modes, for a given lattice volume Vtot- The exact zero 
modes are related to the total topological charge Q top- 
On average, Qtop = 0. The value of Q 2 op fluctuates, 
however. The instanton-like picture presumes that there 
are lumps of topological charge close to Qt op = ±1 oc- 
cupying sub- volumes of order ^-qq D ■ There are no fluc- 
tuations of the topological charge on the scale of the 
lattice spacing a. Indeed, instantons are topological ex- 
citations with smallest action possible and the probabil- 
ity to find instanton of size of order a is proportional to 
a high power of a. Quantitatively, the strength of the 
topological fluctuations is related to the r\ mass 0: 



Xl 



(Q 2 



top/ 



2N 



f 



(213 MeV) 4 



(10) 



where the numerical value is borrowed from Ref |15j . 

There are also near-zero fermionic modes whose total 
number is proportional to the total volume and which 
are associated, in the zero approximation, with the orig- 
inal lumps of the topological charge, see above. Since 
the modes are nearly degenerate, the interaction be- 
tween instantons results in the derealization of these 
zero modes, for details see, e.g., 0. Quantitatively, this 
generic feature of the instantonic vacuum is manifested 
through the Banks-Casher relation [Hj: 



(Vn/>) = -7rp(A„^0) 



(11) 



where p(\ n — > 0) is the density of the (delocalized) zero 
modes in the limit of infinite volume Vtot- For finite 
lattice volumes these modes are near-zero. 

Thus, Eqs I|1U|) and (|llfl predict that the number 
of exact zero modes is proportional to \JVtot while the 
number of near-zero modes is proportional to the total 
volume, Vtot- Moreover, the instantonic picture fixes 
the scaling laws according to which the corresponding 
coefficients of proportionality arc in physical units and 
not dependent on the lattice spacing a. 



were actually many studies of the Witten-Veneziano and 
Banks-Casher relations on the lattice. However they 
used mostly such versions of the lattice fermions that 
the topological charge has intrinsic ultraviolet noise. To 
avoid this, one has to modify, smoothen the original 
gauge fields which fluctuate on the lattice spacing a, for 
review see, e.g., |TT]- As a result, no dependence on 
a could actually be measured. The use of the overlap 
fermions allows to study original gauge fields configu- 
ration, see, in particular, |12j . Study of dependence on 
the lattice spacing a of various observables is one of the 
main objectives of the present paper. 

More explicitly, the massless overlap Dirac operator 
is given by |[J: 

D =-A 1+ vm)- A = D '*- P «- (12) 

where A is the Wilson Dirac operator with negative 
mass term. Anti-periodic (periodic) boundary condi- 
tions in time (space) directions were employed. It 
turns out that the optimal value of p parameter is 
1.4. Furthermore, we have used the minmax polyno- 
mial approximation |13| to compute the sign function 
sign(A) = A/vAAJ = 75sign(_ff), where H — ^A is 
hermitian Wilson Dirac operator. In order to improve 
the accuracy and performance about one hundred low- 
est eigenmodes of H were projected out. Note that the 
eigenvalues of l|12J) lies on the circle of radius p centered 
at (p, 0) in the complex plane. In order to relate them 
with continuous eigenvalues of the Dirac operator the 
circle was stereographically projected onto the imagi- 
nary axis |14| . 

The calculations were performed on two subsets (Ta- 
ble 01 of statistically independent SU(2) quenched con- 
figurations generated with standard Wilson action. For 
the subset A the gauge coupling was chosen is such a 
way that the physical volume remains the same, the cor- 
responding lattice spacings were determined by interpo- 
lating the data of Ref. 15 . The subset B was used to 
determine the IPR volume dependence at fixed spacing 
a = 0.1394(8) fm. 



RESULTS FOR TOPOLOGICAL 
SUSCEPTIBILITY 



4. MEASUREMENTS 



We are using the overlap Dirac operator [§]. The 
advantage of the overlap operator is that it preserves 
the chiral symmetry and allows to study the proper- 
ties of the Dirac modes from first principles JUj . There 



Let us address first the issue of the scaling of topo- 
logical susceptibility \t = (Qtot) /Vtot, where the topo- 
logical charge is defined by the index of overlap Dirac 
operator Qt op = n+ — n_ and n± is the number of ex- 
act zero modes with positive (negative) chirality. The 
topological susceptibility calculated on the subset A 
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a,fm 


L s 


L t 


yphys^ 


N 
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Subset A 


2.3493 


0.1397(15) 


10 


10 


3.8(2) 


300 


2.3772 


0.1284(15) 


10 


14 


3.8(2) 


91 


2.3877 


0.1242(15) 


10 


16 


3.8(2) 


198 


2.4071 


0.1164(15) 


12 


12 


3.8(2) 


179 


2.4180 


0.1120(15) 


12 


14 


3.8(2) 


149 


2.4273 


0.1083(15) 
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16 


3.8(2) 
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2.4500 


0.0996(22) 
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14 


3.8(3) 


200 


2.5000 


0.0854(4) 
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18 


3.92(7) 
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2.3500 


0.1394(8) 


10 


10 


3.8(1) 


100 


2.3500 


0.1394(8) 
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14 


5.3(1) 


100 


2.3500 


0.1394(8) 


12 


12 


7.8(2) 


100 


2.3500 


0.1394(8) 


12 


18 


11.7(2) 


100 


2.3500 


0.1394(8) 


14 


14 


14.5(3) 


94 


2.3500 


0.1394(8) 


14 


16 


16.5(4) 


42 



Simulation parameters. 




0.005 



0.015 



0.01 
a , fm 

Scaling of topological susceptibility with lattice spacing; 
solid curve is the best fit x 1 ^ 4 = Xc/ 4 +ca 2 . Square cor- 
responds to the conventional value (213(3) MeV) 4 , |15| 



The spectral density of eigenmodes of the overlap Dirac 
operator and its interpolation to \ n — > at /3 = 2.5 
on 18 x 16 3 lattice. The value of the condensate (ipip), 
Eq. 1111 . is perfectly consistent with results of Ref. |16|. 




200 250 
X, MeV 



350 400 450 



IPRs for LDEs at various lattice spacings and fixed 
physical volume. The "mobility edge" \ cr « 200 MeV 
is clearly seen. 

modes, see Figure |5| The value of the condensate ('tptp) 
determined through the Banks-Casher relation turns to 
be close to the value obtained earlier with the Wilson 
fermions, see, for instance. |16j. 



scales up to the a 2 corrections (Figure and equals 
to Xo = (225(3) MeV) 4 which is in good agreement with 
the conventional value (213(3) MeV) 4 [T5] . 

The observed scaling of the topological susceptibil- 
ity confirms that the discretization errors and the finite 
volume effects are small. Note that the Wilson gauge 
action is plagued by lattice dislocations and one might 
have speculated that spurious fermionic zero modes are 
generated. However, Figure ^ confirms that the dislo- 
cations are inessential for overlap Dirac operator which 
is insensitive to the ultraviolet noise. 

Moreover, we have measured the density of near-zero 



6. SHRINKING OF THE VOLUME OCCUPIED 
BY FERMIONIC MODES 

A natural measure of the eigenmodes localization 
is provided by the inverse participation ratio (IPR) I\ 
which is defined as follows (for review and applications 
see for example Ref. ^7]). Let p\(x) be the normalized 
bilinear 

x 

where ip\ (x) is the eigenmode of the overlap Dirac oper- 
ator in the given gauge field background with virtuality 
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A, Dip\ = Xip\. Then for any finite volume V the IPR 
I\ is defined by 



h 



(13) 



and characterizes the inverse fraction of sites contribut- 
ing to the support of p\(x). Note that for delocalized 
modes p\(x) = 1/V and hence I\ = 1, while an ex- 
tremely localized mode, p\(x) = S XlXo , is characterized 
by I\ = V. Moreover, for eigenmodes localized on frac- 
tion / of sites (so that sup p\ = Vf = f V) we have 
I\ = V/Vf = l/f. If we allow for a mixture of both 
localized and extended modes the average value of IPR 
is given by: 

(h) =c + Cl V/V f . (14) 

Usually ^7], one considers localization only in terms 
of the total volume. Very recently, there appeared data 
that the localization volume can depend on the lattice 
spacing as well. Namely, such effect was observed first 
in Ref. |18j for fermionic modes (in a particular version 
of lattice fermions, so called Asqtad fermions). Most 
recently localization properties of scalar probe particles 
were investigated ^5] and, again, strong dependence on 
a was observed. Observed dependence on the lattice 
spacing seems specific for dynamics of Yang-Mills theo- 
ries. 

The localization properties of LDEs in our measure- 
ments are illustrated on Figure EH where we plot the in- 
verse participation ratios for modes < A < 450 MeV 
at various spacings and fixed physical volume (subset 
A, Table |2l . One can see that there is a critical value, 
X cr w 200 MeV, above which all the states are in fact de- 
localized and their IPRs are lattice spacing independent. 
However, for small eigenvalues below X cr the value of 
IPR grows and one cannot exclude presence of localized 
states. The mixture of localized and extended modes 
could be quantitatively characterized by the volume de- 
pendence of IPR at fixed lattice spacing, see Eq (|14|l . 
We have computed the average value I\ = J2\(I\) /Ni 
of IPRs for modes < A < A = 50 MeV, where Ni 
is the total number of modes in this interval, and the 
analogous quantity Iq for exact zero modes. Fitting 
the data to Eq. (|14|) we found for exact zero modes 
c = 3.7(5), ci = 0.38(5) with x/ n -d.f — 3, and for low- 
lying modes c = 1.9(2), ci = 2.26(3) with x/n-d.f = 2. 

Moreover, one can see from Figure that the LDE's 
localization degree depends non-trivially on the lattice 
resolution. To quantitatively investigate this depen- 
dence we have plotted Iq and against the lattice spac- 
ing on Figure [SI It turns out that our results are highly 
robust with respect to the actual value of A as long as 
A < 50 MeV so that we could take A = 50 MeV to 




Volume dependence of Io, h at fixed lattice spacing. 




0,11 
a, fm 

Scaling of 7o and 7; with a — > 0. Lines correspond to 
the fitting curves I15t for d = 0, 1, 2, 3. 



improve the statistics. It follows from Figure [S] that the 
inverse participation ratios for both zero and low lying 
eigenmodes seem to be divergent in the limit a — ► 0. 
Note that the average values of IPRs are higher than 
those of Ref. [Hj. 

Furthermore, it is straightforward to estimate the 
dimensionality d of the objects which localize the low 
lying eigenmodes. Indeed, at fixed physical volume we 
have 

h = V/V f = Va 4 / (V/a 4 ) = a d - 4 V phys / Vf ys 
and therefore the lattice spacing dependence is given by 



(h) = b a + b 1 -a 



d-4 



(15) 



We fitted our data to Eq. O for d = 0, 1, 2, 3, the fit- 
ting curves are shown on Figure [SJ As a matter of fact 
the IPR data for low lying modes strongly suggest that 
the dimensionality of underlying objects is zero, d = 0; 
the relevant % 2 /n.d.f. is 2.3 which should be compared 
with its values 3.4 (d = 1) and 4.9 (d = 2). As far as 
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the exact zero modes are concerned our data favor one 
dimensional (d = 1) localization regions; x 2 /n.d./. in 
this case is 0.38 while it is 0.95 for d = and 0.85 for 
d = 2. The error bars are not small enough, however, to 
rule out reliably that the exact zero modes in the limit 
a — * are localized also on point-like objects. 



7. EFFECT OF THE REMOVAL OF THE 
VORTICES 

Thus, our results indicate fine tuning of the 
fcrmionic modes. Namely, their volume shrinks without 
affecting the eigenvalues. We have already mentioned 
(see Eqs. I|6l7[l ) that fine tuning was observed first for 
the central vortices. A natural question arises, whether 
the fine tuning of the fermionic modes is related to 
existence of the fine tuned vortices. A standard way 
to probe such a relation is the so called removal of the 
vortices j23 ED Namely one multiplies the original 
link matrices U^x) by their Zi projected values: 



(16) 



One can show that the modification Hlfijl affects (up to a 
gauge transformation) only a 3d fraction of the total 4d 
volume of the lattice [22] • In other words, this fraction 
is proportional to (a • Aqcd) and vanishes in the limit 
of a — * 0. On the other hand, the vortices are removed. 
For 

The impact of the change Hlfijl on the properties of 
LDEs is illustrated on Figures |HJ Namely, all the 
low-lying modes, with exception of approximately 4% 
of exact zero modes, disappear ^ . For the modes re- 
maining in the spectrum the value of IPR is of order 
unit. In other words, removal of the vortices destroys 
topology-related Dirac modes. 

As is mentioned above the change affects a 3d 
sub-volume of the whole lattice. Further information 
on the manifold crucial for the chirality breaking can 
be provided by observations on the shrinking of the lo- 
calization volume as function of the lattice spacing. If 
indeed, the volume shrinks to points, d=0, as favored 
by the data above, then it is only a small subspace of 
the 3d volume that is related to chiral symmetry. Pure 
geometrically, such points could well be the points of 
self-intersection of the vortices. Theoretically, possible 
connection between self-intersections of the vortices and 
zero modes was considered, within a particular model, 
in Ref. |23) . In principle, this conjecture on the relation 



350 400 



The spectral density of overlap Dirac operator eigen- 
modes for vortex removed configurations. The lattices 
used are 18 x 16 3 at p = 2.5. 




200 250 
X. MeV 



^Note that a similar effect was observed in Ref. 1211 where 
chirally improved lattice Dirac operator was studied. 



IPR for LDEs for vortex removed configurations ob- 
tained at f3 = 2.5 on 18 x 16 3 lattices. 



between self-intersections of the vortices and localiza- 
tion of (near) zero modes could be investigated through 
direct observation. At the moment, however, such data 
are not available. 



8. CONCLUSIONS 

To summarize, we have studied properties of low- 
lying Dirac modes in the vacuum state of SU(2) glu- 
odynamics using overlap fermions. We work with the 
original gauge fields, which fluctuate on the scale of the 
lattice spacing, without any smoothening. The Witten- 
Veneziano and Banks-Casher relations predict scaling 
laws for densities for exact and near zero modes, re- 
spectively. Namely, these densities should depend only 
on the physical scale, Aqcd- From the point of view of 
the lattice measurements, these are strong constraints. 
The results of the measurements do agree with the the- 
oretical expectations. 
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A novel feature uncovered by the lattice simulations 
is that the low-lying modes are localized on the volumes 
which shrink as a power of the lattice spacing a. For ex- 
act zero modes the power is close to three, Vx=o ~ » 3 
while for near-zero modes it is rather four, V\— >o ~ a 4 , 
for error bars see Figure 

Moreover, the removal of the central vortices from 
the original configurations eliminates all the topological 
fermionic modes. This observation favors speculations 
that the lumps of the topological charge are associated 
with self-intersections of the vortices. 

Note that even if this picture is correct, this does 
not automatically imply that the instantonic picture is 
wrong. There could be two dual pictures of the same 
effect, valid for measurements with various resolutions. 
In terms of soft fields, instanton-like fluctuations could 
dominate. In terms of fields measured with fine reso- 
lution, of order a, the fine-tuned surfaces could be re- 
sponsible for the lumps of the topological charge. It goes 
without saying that the hypothesis on the crucial role 
of the surfaces is to be scrutinized much further in the 
lattice measurements. If existence of the dual pictures 
is confirmed by measurements, one could speculate that 
this, empirical duality reflects duality on the fundamen- 
tal level. The vortices themselves could be identified 
with dual strings 0], while their self- intersections could 
be basic element in the mechanism of the chiral symme- 
try breaking in the dual picture. 

One cannot rule out that the lattice-spacing depen- 
dence we found will change on finer lattices. In princi- 
ple, it could happens that the ~ a -4 dependence of IPR 
shown on Figure will eventually flatten in the limit 
a — > so that the LDEs localization will identify some 
4-dimensional objects of finite physical extent. Note 
however, that the characteristic size of these objects is 
definitely less than 0.1 fm and likely to be much smaller 
than that. Although this scenario is not excluded a pri- 
ori and requires further investigations, it seems artificial 
to our mind. 
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